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Abstract 

By using perturbation theory we derive an expression for the elec- 
trical field of a Gaussian laser pulse propagating in a gas medium. 
This expression is used as a trial solution in a variational method to 
get quasianalytical solutions for the width, intensity and self-focusing 
distance of ultrashort pulse. The approximation gives an improved 
agreement with results of numerical simulations for a broad range of 
values of the input power of the pulse than previous analytical results 
available in the literature. 

1 Introduction 

The propagation of high-peak power femtosecond laser pulses in optical me- 
dia has attracted considerable attention recently [HE]- These pulses undergo 
dramatic changes in their temporal, spatial and spectral properties due to 
their nonlinear interaction with the medium. The most fundamental process 
is self-focusing which causes the pulse to be compressed in space, result- 
ing in an increase of the peak intensity Self- focusing would result in 
a catastrophic collapse of the pulse, but the self-focusing process is usually 
balanced by multiphoton ionization or excitation of the atoms or molecules 
in the medium, since the resulting (quasi-)free electron density defocuses the 
pulse. 

The analysis of optical pulse propagation is usually based on the de- 
scription of the pulse in terms of its complex field envelope, neglecting the 
underlying rapid oscillations at its carrier frequency. The resulting slowly 
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varying envelope approximation reduces the Maxwell's equations to higher 
dimensional nonlinear Schodinger equations (NLSE). These equations are 
not integrable, so they do not have soliton solutions. However, they possess 
stationary solutions which are unstable on propagation. The self-focusing 
distance is determined as the (first) point of infinite intensity in the solution 
of the NLSE. 

Most of the quantitative analysis of self-focusing (and further propaga- 
tion) of the pulse results from numerical computation. For example, from a 
curve-fitting procedure based on numerical simulations a popular analytical 
formula for the location of the first singularity (self-focusing distance) for the 
propagation of a CW laser beam has been given by Dawes and Marburger 



where V is the input power scaled in units of the critical power for self- 
focusing, P crit = Ao/27rnon 2 , with Ao is the wavelength and no and ni are 
the linear and nonlinear indices of refraction, respectively [H]. Formula (JI| 
has been verified in many experiments and numerical simulations. Besides 
its accurateness and usefulness for estimations it however does not give much 
insight in the physics of the underlying process. 

In order to get a deeper insight into a physical process, it is often useful 
to use approximative theories, even if full numerical solutions are available. 
A few approximations have been discussed already in the early review pa- 
per on self-focusing by Marburger The resulting expressions provide, 
however, rather qualitative than quantitative estimations only. Alternative 
approaches have been proposed to analyze the effect of self-focusing, namely 
a systematic perturbation theory P E] , a ray-equation approximation [3] , an 
approach based on Fermat's principle [S], a variational analysis j^j, a paraxial 
parabolic approximation [TU], or a source-dependent expansion method [TT] . 
The predictions of these models for the self-focusing distance agree quanti- 
tatively well with the exact solutions from numerical simulations for input 
powers close to the critical power, but deviate quickly as the power increases. 

In this paper we present a quasianalytical approximation which provides 
a description of the self-focusing phenomenon not only qualitatively but also 
quantitatively. It is based on the variational approach using a trial solution, 
which contains a first-order perturbation correction of the phase resulting 
from the transverse distortion of the pulse. As we show below, predictions 
for the self-focusing distance, the on-axis peak intensity and the radius of 
the pulse within this ansatz are in better agreement with those of numeri- 
cal simulations over a broad range of input powers than previous analytical 
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results available in the literature. We show how the pulse separates into an 
inner and an outer self-trapped component along the propagation distance. 
The inner component self-focus, while the outer one stays as background. 

The paper is organized as follows: First, we derive the phase correction 
of the pulse using first-order perturbation theory. The corrected form of 
the pulse is then used as trial solution in the Lagrangian of the system to 
obtain equations for the self-focusing distance, the on-axis peak intensity 
and the radius of the pulse. Next, the predictions of the approximation will 
be compared with those of the earlier models and of numerical simulations. 
Finally, an analogy with the problem of a particle in a finite quantum well 
is presented. In this analogy, the bound state and its exponential-decaying 
wings penetrating a smooth potential wall correspond to the inner and outer 
component of the pulse, respectively. The self-focusing is associated with the 
shrinking of the width of the quantum well. 



2 Variational approach 

We study the propagation of a linearly polarized laser beam in gases. For the 
explicit calculations we do not restrict to any specific gas. Our analysis is 
based on the scalar wave equation, which can be obtained from the Maxwell 
equations (e.g. fHj): 

dlE + AE- = ^^d?E 3 + ^d t J. (2) 

Here, x is the third order nonlinear susceptibility coefficient of the medium, 
eo is the linear dielectric constant of the gas and c is the speed of light. 
Since we are interested in the dynamics of the pulse up to the self-focusing 
distance only, we neglect the plasma effects (J = 0), which usually balance 
the catastrophic collapse and, hence, influence mainly the propagation of 
the pulse after the self- focusing point. Using the slowly varying envelope 
approximation, Eq. (j2j) can be written in the retarded coordinate frame 
(t — > t — z/vg) as: 

■a i 1 A | | 1 2 

lOcU-\ — Aj_ u + - — \u\ u 

-i£(*" + i§*») = °- < 3 > 

where u = £/v^0) £ = z/L D , r = t/T , and the transverse coordinate r 
is given in units of the length of the pulse w (w is the radius at 1/e 2 of 
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Figure 1: Results of numerical simulations for the normalized pulse shape 
for P/P cr = 5: f = (dotted-dashed line), f = 0.15 (dashed line), f = 0.20 
(dotted line), f = 0.23 (solid line). 



the irradiance). The characteristic lengths are given by Lp = k z w\j2^ Ld = 
T 2 /2/c 2 , L d > = T 3 /2fc 3 and L nl = (n 2 &o^o) -1 , where k Q = 2tc/X, k z = n k . k 2 
and &3 are the second and third order group-velocity dispersion coefficients, 
respectively, T is the duration of the pulse (T FWHM = ^2 log2T ), and 
Iq is the input peak intensity (in units of W/cm 2 ). We further neglect in 
Eq. Q the effect of the group velocity dispersion to get the well-known 
(2+l)-dimensional NLSE, 

id^u + ^A ± u + F[uu*]u = 0, (4) 

where the star stands for complex conjugation and F[uu*] = ci\u\ 2 with 
c\ = Lo/Lni. During the early stage of propagation of the pulse (up to 
the self-focusing point) in gases the effect of the group velocity is negligible, 
namely if L n i <C < L$. Thus, for the present analysis we expect that the 
condition is fulfilled as long as the self-focusing distances are small, i.e. as 
larger the input power of the beam is. 



2.1 Phase correction 

The variational method has been applied to the (2+l)-dimensional NLSE 
recently [H] to study qualitatively the phenomenon of self-focusing in air. As 
mentioned at the outset the method is as good as the trial solution. Usually 
the self-similarity assumption is used in order to integrate over the transversal 
coordinates taking into account that the total power of the NLSE is invariant. 
This is done by assuming that the shape of the pulse remains unchanged up 
to the self-focusing point. The assumption preserves the lens transformation 
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of the NLSE. However, from numerical simulations one can observe that the 
initial pulse shape gets distorted during the self- focusing process (c.f. Fig. [TJ) . 
In fact, the pulse separates into two components as it propagates [Zj, 



where i/jq is the high intensity inner core of the pulse, which self-focuses and 
eipi is the low intensity outer part, which propagates forward following the 
linear propagation mode. In recent experiments and numerical simulations 
[T^l FEU ITU Ho"] it has been shown that this weak large background plays a 
decisive role in the propagation and filamentation process. 

We expect that the self-focusing process is influenced by the interaction 
between the inner core and the outer part too. To this end we use a trial 
solution in the variational method, in which the perturbation by the back- 
ground is taken into account as a correction in the phase. We first separate 
amplitude and phase of the pulse as: 



(5) 



\u{£,r)\exp(iS(£,r)) 
|V"b(f,r)|exp(i5b(f,r)) 
+<#i (£,?") | exp (iSi(£,r)) 



(6) 



and by substituting in Eq. 



(JH) we get up to order e that 
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Here, 



no&r) = F[|^,r)| a ]--(^,r)) 2 



(9) 
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and 



m (£, r) = e| f 2F'[|^ (e, r) | 2 ] |Vb(£, r) 1 r) | 
|^(e,r)|A ± |^ (^r)| + A ± |^(e,r)| 



4|^ (e,r)| 2 4|Vo(£,r)| 



4|^ (£,r)| 
x cos (5'o(^r)-S' 1 (e,r)) 

( \MM\a ± Si{M 

V 4|^ (e,r)| 
^i(g,r)a r 5i(e,r) 
2|Vo(£,r)| 

x sin(5 (e,r)-5i(e,r))^ (10) 



The desired trial solution is then given by 

= ^(£,r)exp(iSi(f,r)), (11) 

with Si(0, r) = 0. Here the amplitude corrections of order e have been 
neglected, while phase corrections have been kept. Note that «(£, r) still 
fulfils the self-similarity assumption, namely |w(£,r)| 2 = iV'o(^) r )| 2 - 

The phase correction Si(£,r) is obtained using first-order perturbation 
theory as follows. Inserting Eq. (jSj) into Eq. (J3J and collecting powers of e 
we get 

e° : ^o + ^A ± ^ + F[^o]^o = 0, (12) 

e 1 : z%V>i + ^A ± V>i + ^[VW>o]V>i 

+F> ^oWi + i^oWoVi = 0, (13) 

where the prime stands for the derivative with respect to the argument. Note, 
that the solution of Eq. (|13J) has the form 



ipi = d a (t)ipo, (14) 

where a(£) is the length of the pulse. Thus, the phase correction can be ob- 
tained for any shape of the inner core by taking its derivative and separating 
the phase. 
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We do not consider the presence of any external perturbation, like losses or 
ionization. However, the theory presented here can be extended considering 
these terms following the same lines as in |J] . We have checked carefully that 
the predictions of the numerical calculations for the self-focusing distance do 
vary by less than 1%, if we include external perturbations. Other properties, 
such as the pulse length or the on-axis intensity (up to the self-focusing point) 
are even less sensitive to these effects. It is therefore justified, for the sake of 
simplicity, to neglect external perturbations. 



2.2 The Gaussian pulse shape 

As outlined above, the phase correction can be obtained for any shape of 
the pulse. In order to investigate the effect we apply the theory to the most 
important case of a Gaussian pulse shape below. In this case 

where p = r/a(£) and bo is a constant describing the initial wave front di- 
vergence of the pulse (i.e. b = for collimated pulse, or b = —1/f for the 
case of an external lens with focal length /), we obtain the phase correction 
from the solution of Eq. (JIT ~ 



^l(6; r ) = ^77T ex P(-p 2 ) ex P(* & oa(Op 2 ) 



x\2p i -ib ,a(t)p i -lj, (16) 
as follows. Substituting Eqs. (fTKj) and (|TK|) into and (fTTljl we get 

MM = -\(d r S(Z,r)) 2 



+ ^7T(p 2 -l + A 2 ciexp(-2 /3 2 )) 



(17) 



and 



n ±(M = (2p 2 exp(-2p 2 ) -exp(-2p 2 )) 
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where a(£) is an undetermined function. In Ref. jH] it has been shown that 
a(£) can be determined using 

/5 = -a 3 (Oa"(0, (19) 

which is proportional to the excess power, as long as the excess power is 
small. With this relation one is able to solve the integral (JHJ) in terms of (3 for 
Gaussian pulses. However, the solution of the integral is ill-posed, since the 
integrals in the solution diverge. Therefore, we have adopted an adiabatic 
approximation instead, namely that the transverse form of r) remains 
unaffected along the propagation distance, i.e 

Si((,r) = ni(£,r)£ (20) 
2 A 2 Cl e«e 



2ee 



(2p 2 exp (-2p 2 ) - exp (-2p 2 )) 



(l-2p 2 ). (21) 



The last term of Eq. (|21jl arises due to the diffraction term in Eq. (fT3*|) for 
the outer part of the pulse. This term should not effect the propagation of 
the inner core part and we neglect it. The validity of this assumption can be 
shown using the variational method (see Appendix El). Thus, finally we get 

u(£,r) = -J— Aexp(p 2 ) 

x exp(i& a(Op 2 )exp(i,Si(f,r)) (22) 

with 

Si(^r) = 2A l%** (2f - l)exp(-2p 2 ). (23) 

2.3 Lagrangian and variational analysis 

The Lagrangian functional for the system (J1J) reads, 

£ = ^{ud^u* -u*d(:u) + ^|%x| 2 - ^ci|m| 2 . (24) 



From Eqs. f!22[) and (|46|l we define 



u(f,r) = — |yA(0exp (-p 2 ) exp («6 a(Op 2 ) 

x exp (20(0) exp (iSi(£,r)) (25) 
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with 



Si{M = c(Oa(0 (2p 2 - 1) exp (-2p 2 ), (26) 

where the variational parameter A(£), c(£) and </>(£) are sufficient to 

describe the dynamics of the problem. Here, we consider a collimated input 
laser pulse, focused by a lens, given by 

u(0, r) = uo exp (— r 2 ) exp (— r 2 ) exp f — i~^J ■ (27) 

where / is the focal length of the lens in units of diffraction length Lp. 
Note that the temporal dependence of the variational parameters in (J23j) 
is not written explicitly. Inserting Eq. into the Lagrangian (j2H and 

integrating over the transverse coordinate we obtain the reduced Lagrangian 



(£) = 2tt J Crdr, (28) 
o 

which depends on the variational parameters and the independent variables 
£ and r only. The equations of motion for the variational parameters are 
given by 

^d {d ^ ) (C)-dM = 0, (29) 

where \ii= A, a, c, 4> (i=l,...,4). This leads to the following set of coupled 
equations: 

a'(£) = y c(£) + &q, (30) 



V(0 = 0, (32) 

oo 

where V(g) = P(0/ p cr with P(£) = 2vr /|u|rdr = ttA 2 (£)/2, and the 

o 

critical power P cr = tt/ci. 

From Eqs. pH|l - ([3*2^ with the initial conditions, a(0) = 1 and c(0) = 0, 
we get 

32 

a 3 (£K(fl = -0 = --g(P-l)- (33) 
9 



or 



i(0 = 6o0 2 -/?e (34) 

(35) 



The self-focusing distance, 

1 

is obtained for the collapse of the pulse (a(£) = 0). Note that 



£./ = -^s— r, (36) 



+ i (37) 



e./ C/(/'n = 0) / 

in agreement with the lens transformation property of the NLSE The 
on-axis intensity is given by 

,2 ^(0 



m = mm- a2(e) 

^ (3. 



(1_ 6o0 2 -^ 2 

3 Results and comparisons 

Below we present the predictions of the present theoretical ansatz and com- 
pare them with those of other semi-analytical estimations and the results of 
numerical calculations using Eq. (J3J . The comparisons are performed for the 
case of a collimated Gaussian pulse without external focusing. We should 
note that this is the most extreme case. The presence of a finite external 
focusing reduces the relative errors between the theory and simulations. 

First, we concentrate on the results of the present ansatz for the self- 
focusing distance, £ s / (c.f. Eq. (|36|) ). In Fig. |2]the relative error (solid line), 

t - £ ( M ) 

Vsf = (39) 
W 

with respect to Marburger's formula, Eq. (JTj), is shown as a function of the 
input power (scaled in units of the critical power). Note that Eq. has 
been derived from numerical simulations and is, hence, an excellent estimate 
of the exact self-focusing distance. The comparison shows that the relative 
error is largest (about 40%) for input powers close to the critical power and 
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Figure 2: Relative error T] s t of the predictions of the present theory for the 
self-focusing distance (solid line) with respect to the Marburger's Formula as 
function of V = P/P cr . Also shown are the relative errors for earlier models 
(see text). 



is less than 20% for 2.5 < V < 20. The large error near the critical power for 
self-focusing might be due to the neglect of the group velocity in the analysis. 

Also shown in Fig. El are the relative errors (with respect to Marburger's 
formula) for the predictions of the self-focusing distance resulting from analy- 
ses by means of either one of the following methods, ray-equation approxima- 
tion [3], Fermat's principle jH], variational analysis j^j, paraxial parabolic ap- 
proximation [TU], or source- dependent expansion method [TT] (dashed line): 

and by using perturbation theory (dotted line, [Zj): 

(41) 




N C VV^T 

where the constants M = 0.55 and iV c = 1.86 are derived from the Townes 
soliton shape. The value of ^) in Eq. (}4~T|) is given in diffraction units and 
therefore four times larger than in [7j. The latter formula, Eq. 1)41)1. is valid 
for V < 2. 

Both estimations are in excellent agreement with the results of Mar- 
burger's formula near the critical power but strongly deviate as the power 
increases at larger powers. It is clearly seen, that already at V ~ 2 the predic- 
tions within the present theory are in closer quantitative agreement with the 
exact results. For example, at V — 10 the predictions of the earlier theories 
deviate from the exact result by more than 70% in the case of perturbation 
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Figure 3: Comparison of the predictions of the present theory (solid lines) 
for the pulse length with numerical simulations (dashed lines) for V = 3 (A), 
V = 5 (B), V = 10 (C), V = 20 (D). 



theory (see Eq. (|41jl ) or by more than 100% in the case of other theories 
(see Eq. (|4"UJ) ). while the present theoretical prediction agree within about 
20% with Marburger's formula. We attribute this to the phase correction 
obtained from the interaction with the weak large background, which has 
not been taken into account in the earlier work. In fact, from our ansatz we 
observe that the background does not diffract and receives energy from the 
inner part of the pulse during the self-focusing process, as we will discuss 
below in Section 0] We also mention that the expression in Eq. (|4"Uj) for 
the self-focusing can be obtained by taking into account the two first leading 
terms of the power series of the phase correction in Eq. (j2*Uj) only. This 
means that the earlier ansatz (see for instance jHJ) is a limiting case of Eqs. 
1)25) and (J2EJ). 

We find also that the estimations of the present theory and the numerical 
simulations for the pulse length and the on- axis intensity agree within 30%. 
In Fig. El a comparison of the pulse length a(£), Eq. (J34"|) (solid line), and the 
l/e 2 -level width of the intensity profile using numerical simulation (dashed 
line) of a collimated Gaussian pulse for different input powers is shown. 

Since the range of values of the ordinate in Fig. El is the same (but not 
those of the abscissa) we have evaluated an relative error using the inverse 
functions, 

£ [^theory] £[Osjmu/<jtion] (42) 

C [Q'simulation\ 

The results are presented in Fig. EJa) as a function of a(£) for several values 
of the input power. We observe that rji is in general lower than 10% for 
a(£) ~ 1, i.e. in the early stages of the propagation (c.f Fig. El). Close to 
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Figure 4: Relative errors 771 (panel a) and 772 (panel b) for different input 
powers, V = 3 (solid line), V = 5 (dotted line), V = 10 (dashed line), 
P = 20 (dotted-dashed line). 
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Figure 5: Comparison between the intensity profile of the pulse (solid line) 
obtained from the simulation and the trapping potential derived from the 
theory at £ = 0.2 with V = 5. 



the self-focusing point, a(£) ~ 0, rji increases to above 20% for high input 
powers (V > 10). This shows that in the simulations the pulse undergoes a 
stronger deformation than predicted by the present ansatz. 

In the same way we define the relative error for the on-axis intensity, 
namely 

£ [•^tft.eorj/J 'simulation] l'AQ\ 

V2 = 7?y ; • (4dj 

S l 1 simulation] 

As can be seen from the results in Fig. EJb) 772 is lower than 20% for V < 10. 
In general, 772 tends to be constant for larger values of the intensity. 



4 Quantum well analogy 

In order to further understand the role of the weak large background on the 
dynamics of the pulse, we have investigated our ansatz using the quantum 
well analogy. In this picture £ acts as a fictitious time variable, w(£,r) as a 
bound state and —</>'(£) as the energy function of the system. By substitut- 
ing Eq. ()25|) into Eq. (jlj) and taking the real part, one finds an equation for 
—</>'(£). If we neglect the diffraction part of this equation (neglecting kinetic 
energy in the quantum analogue), it is straight forward to derive a potential- 
like function. In Fig. EJthe normalized shape of this potential is compared 
with the normalized shape of the intensity profile of the pulse obtained from 
numerical simulations. We observe that the shape of the potential is a well 
which contains two regions, an inner and an outer part. The inner part ex- 
tends from the center of the well (r = 0) up to the absolute minimum. The 
outer part includes the rest, namely the potential wall which smoothly grows 
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Figure 6: The upper panel contains an example of the intensity distribution 
obtained from simulation (solid line) for V = 5 and £ = 0.2. The distribution 
can be fitted by the superposition of two Gaussians (dashed line). One for 
the inner part, and the other for the background (dotted line). The lower 
panel contains the result of the operation in Eq. (j44j) . Note that the region, 
where the dashed line tends towards the dotted line (inner and outer part 
join), is located where the maximum of the function r\ is. 

radially up to a certain constant value. Notice that the wall is infinitely thick, 
so no analogue of the tunnel effect is possible here. These features of the po- 
tential well reveal already that the pulse possesses two components. In fact, 
we observe in Fig. El that the inner component of the pulse is bell-shaped, 
while the outer component is an analogue of the exponential-decaying wings 
of the probability density of a bound state penetrating a smooth potential 
well. The potential well shrinks radially along the propagation distance and 
the local maximum in the center decreases and eventually disappears at the 
self-focusing distance. So, the inner component of the pulse tends to concen- 
trate at the origin and eventually self-focuses. 

We have estimated the region in which the inner part of the pulse is con- 
centrated as follows. From the intensity profile of the pulse, obtained from 
numerical simulations, one can notice that it is well fitted by the superpo- 
sition of two Gaussian functions, as shown in the upper panel of Fig. |U1 In 
order to determine quantitatively where the inner and outer part join, we 
have performed the following operation in our simulations: 



where |ii(£,r)| 2 is the intensity distribution. Note that rj is constant when 




(44) 
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Figure 7: Comparison of the position of the minimum of the potential (dashed 
line) and the position of the deformation of the Gaussian pulse (solid line) 
vs. £ for V = 5. 
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Figure 8: Estimated fraction of energy, fi, that self-focuses for V — 5. Solid 
line: simulation, dashed line: theory. 

ii(£, r) is a pure Gaussian function. For the present case (c.f. upper panel of 
Fig. r\ has its absolute maximum in the region where the inner and outer 
part of the pulse join (c.f. lower panel of Fig. EJ). Thus, we determine the 
boundary position, rj oin , between the inner and outer part of the pulse by 
the absolute maximum of rj. 

In Fig. [7|we show the comparison between rj oin from simulations (solid 
line) and the position of the absolute minimum of the potential obtained 
from the quantum well potential (dashed line c.f. Fig. EJ) Note that at £ = 
no deformation is observed in the simulations since the initial pulse shape is 
Gaussian. For £ > a deformation of the pulse appears at the origin (r = 0) 
and moves quickly to r ~ 0.65. Afterwards rj oin shrinks and disappears at 
the self-focusing distance. This qualitative behavior is found for all the cases 
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considered in this paper. Thus, our theoretical model does not predict the 
appearance of the deformation at the origin, but the position of the absolute 
minimum of the potential well follows closely the behavior of Tj i n over the 
main part of the propagation distance. 

From the radii shown in Fig0 we can calculate the fraction of energy \i 
stored in the inner component of the pulse, which is compared in Fig. [HI with 
those from the numerical simulations. We observe that there is a energy flow 
from the inner part to the outer part when the pulse shrinks. Note that 
the inner component of the pulse releases more energy in the second half 
of its propagation towards the self-focusing distance, where the shrinkage of 
Tj oin is more pronounced. At the self-focusing distance most of the energy is 
stored in the outer component of the pulse which stays with the pulse and 
does not diffract. We also note that for low powers (V < 5) the present 
model overestimates the release of energy towards the outer region, while for 
higher powers (V > 5) it is underestimated. This difference is related to the 
discrepancy found for the estimation of the self-focusing distance (see Fig. 
EI)- We note, that the ansatz in earlier works (see for instance |9J) would 
lead to a potential well too, however the wall of this potential well has a 
finite thickness. Thus, in the previous ansatz the outer part of the pulse 
diffracts out in contrast to the present model where the pulse is completely 
self-trapped. 

5 Conclusion 

In conclusion, we have revisited the problem of self-focusing of a laser pulse in 
a gas Kerr medium modeled by the (2+l)-dimensional nonlinear Schrodinger 
equation. We have defined a new trial solution taking into account that 
during the self-focusing process the pulse splits in an inner and an outer 
component. In the trial solution the outer component is taken into account 
via a phase correction. We have used a quantum well analogy to explain 
the dynamics of the laser pulse. According to this picture, the laser pulse is 
represented by a bound state of a particle trapped in a quantum well. The 
inner component of the laser pulse corresponds to that part of the bound state 
where the energy of the system is larger than the bottom of the potential well. 
The outer part is the analogue of the exponential-decaying wings of the bound 
state penetrating a smooth potential wall. Here, the self-focusing process of 
the laser pulse is associated to the shrinkage of the width of the quantum 
well. Theory and simulations show that during the self-focusing process the 
energy of the inner component flows to the outer one which stays with the 
pulse and does not diffract. Finally, only a fraction of the pulse energy 
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self-focuses. The present theory provides better agreement with results of 
numerical simulations than most of the theories used in the last decades for 
the self-focusing process of laser pulses. Comparison of the prediction of the 
present theory with Marburger's formula for the self-focusing distance shows 
an agreement within an error of 20% over a broad range of input powers of 
the pulse. Further comparisons with numerical simulations show that the 
pulse length as well as the on-axis intensity are predicted correctly within a 
20% error too. We may emphasize that the comparisons are performed for 
collimated laser pulses, which is the worst case possible. Consideration of 
an external focusing would give even lower errors. Finally, we note that the 
present theory can be easily extended to consider losses and ionization of the 
gas medium. 

A Analysis of the diffraction term 

From Eqs. ([lljh (|15|) and (|21|) we can write 

r) = -^jMO ex P (p 2 ) ex P (* b o a (0 p 2 ) 

x exp ^c(Oa(0(2p 2 exp (-2p 2 ) - exp {-2p 2 ))\ 

x ex p(^( 1 - 2 ^))' ( 45 ) 

where A(£), a(£) and c(£) are assumed to be variational parameters, and k 
is an unknown constant. By doing variational approximation similar to that 
one in sectional but using the ansatz (|45|). one can show that in the critical 
regime (P = P cr ) 

k = constant x a"(0) = 0. (46) 

Thus, in principle we can neglect corrections on the diffraction term and keep 
only the corrections in the nonlinear term in ([45)1 . 
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